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We determine the fraction of tf events with spin correlation, assuming that the spin of the top 
quark is either correlated with the spin of the anti-top quark as predicted by the standard model 
or is uncorrelated. For the first time we use a matrix-element-based approach to study tt spin 
correlation. We use tt — > W + b W~b — > £ + i>bl~i'b final states produced in pp collisions at a center 
of mass energy yfs = 1.96 TeV, where I denotes an electron or a muon. The data correspond to an 
integrated luminosity of 5.4 fb~ and were collected with the DO detector at the Fermilab Tevatron 
collider. The result agrees with the standard model prediction. We exclude the hypothesis that the 
spins of the tt are uncorrelated at the 97.7% C.L. 

PACS numbers: 14.65.Ha, 12.38.Qk, 13.85.Qk 



While top and anti-top quarks are unpolarized in tt 
production at hadron colliders and their spins cannot 
be measured directly, their spins are correlated and this 
correlation can be investigated experimentally The 
standard model (SM) of particle physics predicts that top 
quarks decay before fragmentation [2[ , which is in agree- 
ment with the measured lifetime of the top quark Q. 
The information on the spin orientation of top quarks 
is transferred through weak interaction to the angular 
distributions of the decay products 0, @ • 
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We present a test of the hypothesis that the correla- 
tion of the spin of t and t quarks is as expected in the SM 
as opposed to the hypothesis that they are uncorrelated. 
The spins could become decorrelated if the spins of the 
top quarks flip before they decay or if the polarization 
information is not propagated to all the final state prod- 
ucts. This could occur if the top quark decayed into a 
scalar charged Higgs boson and a b quark (t — > H + b) 0- 
i- 

Recently, the CDF Collaboration has presented a mea- 
surement of the tt spin correlation parameter C in 
semileptonic final states from a differential angular dis- 
tribution @. The spin correlation strength C is defined 
by d 2 a / d cos 6id cos 6% = er(l — Ccos#i cos02)/4, where 
a denotes the cross section, and d\ and O2 are the angles 
between the direction of flight of the decay leptons (for 
leptonically decaying W bosons) or jets (for hadronically 
decaying W bosons) in the parent t and i rest frames and 
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the spin quantization axis. The value C = +1 (—1) gives 
fully correlated (anticorrelated) spins and C — corre- 
sponds to no spin correlation, while the NLO SM predic- 
tion using the beam momentum vector as spin quantiza- 
tion axis is C = 0.777±g;g|| 4]. The DO Collaboration 
has performed two measurements of C in dilepton final 
states [13, EH , where the second analysis uses the same 
dataset as this measurement. None of the previous anal- 
yses has sufficient sensitivity to distinguish between a hy- 
pothesis of no correlation and of correlation as predicted 
by the SM. 

In this Letter, we present the first measurement of spin 
correlation in tt production using a matrix-element-based 
approach, exploring the full matrix elements (ME) in 
leading order (LO) Quantum Chromodynamics (QCD). 
We extract the fraction / of tt candidate events where the 
tt spin correlation is as predicted by the SM over the total 
number of tt candidate events assuming that they consist 
of events with SM spin correlation and of events without 
spin correlation. We use tt event candidates with two 
charged leptons in the final state, where the charged lep- 
tons correspond to either electrons or muons, in a dataset 
of 5.4 fb _1 of integrated luminosity that has been col- 
lected with the DO detector at the Fermilab Tevatron pp 
collider. With a matrix-element-based approach, we use 
the full kinematics of the final state to improve the sensi- 
tivity with respect to using only a single distribution by 
almost 30%. 

The DO detector [l2| comprises a tracking system, a 
calorimeter, and a muon spectrometer. The tracking sys- 
tem consists of a silicon microstrip tracker and a central 
fiber tracker, both located inside a 2 T superconducting 
solenoid. The system provides efficient charged-particle 
tracking in the pseudorapidity region |?7dct| < 3 [13j . The 
calorimeter has a central section covering \rjdct\ < 1-1 
and two end calorimeters (EC) extending coverage to 
|*7det | ~ 4.2 for jets. The muon system surrounds the 
calorimeter and consists of three layers of tracking detec- 
tors and scintillators covering r/det | < 2 [lij , A 1.8 T 
toroidal iron magnet is located outside the innermost 
layer of the muon detector. The integrated luminosity 
is calculated from the rate of inelastic pp collisions, mea- 
sured with plastic scintillator arrays that are located in 
front of the EC. 

We use the same selection of It (ee, e/z, and /z/z) events 
as described in Ref. fill ], therefore only a short overview 
of the selection is given. To enrich the data sample in 
tt events, we require two isolated, oppositely charged 
leptons with pt > 15 GeV and at least two jets with 
Pt > 20 GeV and |?7dct| < 2.5. Electrons in the cen- 
tral (|%et| < 1-1) and forward (1.5 < r/dct| < 2.5) re- 
gion are accepted, while muons must satisfy |?7det| < 
2. Jets are reconstructed with a mid-point cone algo- 
rithm [l5[ with radius 1Z = 0.5. Jet energies are cor- 
rected for calorimeter response, additional energy from 
noise, pileup, and multiple pp interactions in the same 
bunch crossing, and out-of-cone shower development in 
the calorimeter. We require three or more tracks origi- 



nating from the selected pp interaction vertex within each 
jet cone. The high instantaneous luminosity achieved by 
the Tevatron leads to a significant background contribu- 
tion from additional pp collisions within the same bunch 
crossing. The track requirement removes jets from such 
additional collisions and is only necessary for data taken 
after the initial 1 fb _1 . The missing transverse energy 
($t) is defined by the magnitude of the negative vec- 
tor sum of all transverse energies measured in calorime- 
ter cells, corrected for the transverse energy of isolated 
muons and for the different response to electrons and jets. 
A more detailed description of objects reconstruction can 
be found in [l6j . 

The final selection in the e/z channel requires that the 
scalar sum of the leading lepton pt and the px of the two 
most energetic jets be greater than 110 GeV. To reject 
background in ee and /z/z events, where JpT arises from 
mismeasurement, we compute a Ifir significance which 
takes into account the resolution of the lepton and jet 
measurements. We require the significance to exceed five 
standard deviations. In the /z/z channel, events are fur- 
thermore required to have Jf^T > 40 GeV. 

The tt signal is modeled using the mc@nlo [l7j event 
generator together with the CTEQ6M1 parton distribu- 
tion function (PDF) [l8j], assuming a top quark mass 
rat = 172.5 GeV. We generate tt Monte Carlo (MC) 
samples with and without the expected spin correlation, 
as both options are available in MC@NLO. The events 
are processed through herwig [19( to simulate fragmen- 
tation, hadronization and decays of short-lived particles 
and through a full detector simulation using GEANT [2(| . 
We overlay data events from a random bunch crossing 
to model the effects of detector noise and additional pp 
interactions to the MC events. The same reconstruction 
programs are used to process the data and MC simulated 
events. 

Sources of background arise from the production of 
electroweak bosons that decay into charged leptons. 
In the ee, e/z, and /z/z channels, the dominant back- 
grounds are Drell-Yan processes, namely Z/j* — > e + e _ , 
Z/-f* -> t+t~ -> vl + vvl~v, with £ ± = e ± or /z*, and 
Z/7* — > In addition, diboson production (WW, 

WZ and ZZ) contributes when the bosons decay to two 
charged leptons. We model the Z/-f* background with 
ALPGEN [2lJ , interfaced with pythia Q , while diboson 
production is simulated using pythia only. The Z/j* 
and diboson processes are generated at LO and are nor- 
malized to the next-to-next-to-leading order (NNLO) in- 
clusive cross section for Z/j* events and to the next- 
to-leading order (NLO) inclusive cross sections for di- 
boson events [23l l24ll . For all background processes the 
CTEQ6L1 PDF are used. 

Detector-related backgrounds can be attributed to jets 
mimicking electrons, muons from semileptonic decays of 
b quarks, in-flight decays of pions or kaons in a jet, and 
misreconstructed Ifr. These backgrounds are modeled 
with data. Background from electrons that arise from 
jets comprising an energetic ir° or rj particle and an over- 
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TABLE 1: Yields of selected events. The number of tt events 
is calculated using the measured cross section of a t i = 8.3 pb 
and the measured / = 0.74. Uncertainties include statistical 
and systematic contributions. 

tt Z/'y* Diboson Instrumental Total Observed 
341 ± 30 93 ± 15 19 ± 3 28 ± 5 481 ± 39 485 



lapping track is estimated from the distribution of an 
electron-likelihood discriminant in data |16| . In the e/i 
and channels, muons produced in jets that fail to be 
reconstructed can appear isolated. Table Q] summarizes 
the yields for the signal and background contributions. 

To distinguish the hypothesis H of correlated top 
quark spins as predicted by the SM (H = c) from the 
hypothesis of uncorrclatcd top quark spins (H — u) , we 
calculate a discriminant R [25[ defined as 



I? 



P sgD (H = c) 



P ssn (H = u) + P sgn (H = c) ' 



(1) 



where we calculate per-event probability densities, P sgn , 
for tt signal events for both hypotheses constructed from 
the LO MEs M(y,H) [H, 



P sg n(x;H) 



1 



Cobs 



/pdf(<?i) /pDF(<72)d<7idg 2 



(2tt) 4 \M(y,H)\ 2 



W(x,y)d$ 6 - (2) 



Here, cr b s denotes the leading order cross section includ- 
ing selection efficiency, q\ and qi the energy fraction of 
the incoming quarks from the proton and antiproton, re- 
spectively, /pdf the parton distribution function, s the 
center-of-mass energy squared and d$g the infinitesimal 
volume element of the 6-body phase space. The detector 
resolution is taken into account through a transfer func- 
tion W(x,y) that describes the probability of a partonic 
final state y to be measured as x = (pi, . . . ,p n ), where 
pi denotes the measured four-momenta of the final state 
particles. For hypothesis H = c we use the ME for the 
full process qq — >• tt — > W + bW~b — > t + vibt'~v^b aver- 
aged over the initial quarks' color and spin and summed 
over the final colors and spins [26]. For hypothesis H = u, 
we use the ME of the same process neglecting the spin 
correlation between production and decay [261 ] . The tt 
production cross section, a t t, does not depend on the hy- 
pothesis H = c or H — u, and is taken as identical for 
both hypotheses. It is assumed that momentum direc- 
tions for jets and charged leptons and the electron energy 
are well measured, leading to a reduction of the number of 
integration dimensions. Furthermore, the known masses 
of the final state particles are used as input, and it is 
assumed that the tt system has no transverse momen- 
tum resulting in a six dimensional phase space integra- 
tion. More details of the calculation of P sgn can be found 
in 27]. Figure Q] shows the discriminant R for generated 
partons for H = c and H = u for tt MC events. 
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FIG. 1: Comparison of the discriminant R between SM spin 
correlation H — c and no spin correlation H = u at parton 
level. The first and last bin include also the contributions 
from R < 0.29 and R > 0.63. 



To measure the fraction / mea s of events with SM spin 
correlation, we build templates of R distributions for sig- 
nal MC with and without spin correlation as well as for 
each source of background. The templates are compared 
to the R distribution in data and the fraction of events 
with SM spin correlation is extracted. 

In Fig. [5J the measured discriminant R in data is com- 
pared to templates for tt production with SM spin corre- 
lation and without spin correlation including background 
for all dilepton channels combined. The separation be- 
tween H — c and H — u is decreased compared to the 
parton level. 
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FIG. 2: (Color online) The predicted discriminant distribu- 
tion R for the combined dilepton event sample for the fitted 
Oti and /meas compared to the data. The prediction with spin 
correlation (/ = 1) and without spin correlation (/ = 0) is 
shown including background. The first and last bin include 
also the contributions from R < 0.29 and R > 0.63. 



We perform a binned maximum likelihood fit to the R 
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distribution to extract / moas by fitting 

m W = /meas m« + ( 1 - / mC as) m« + £ mf , (3) 
(i) 

where m c is the predicted number of events in bin i for 
the signal template including SM spin correlation, m£ is 
the predicted number of events in bin i for the template 
without spin correlation and mj 1 ' 1 is the sum over all 
background contributions j in bin i. To remove the de- 
pendence on the absolute normalization, we calculate the 
predicted number of events, rrS l \ as a function of f mcas 
and a t i and extract both simultaneously. 
The likelihood function 

N K 

C = Y[P(n®,m®) x Jj0(i*;O,SD fc ), (4) 

i k=l 

is maximized with V(n,m) representing the Poisson 
probability to observe n events when m events are ex- 
pected. The first product runs over all bins i of the tem- 
plates in all channels. Systematic uncertainties are taken 
into account by parameters Vk, where each independent 
source of systematic uncertainty k is modeled as a Gaus- 
sian probability density function, Q (v; 0, SD), with zero 
mean and an rms corresponding to one standard devia- 
tion (SD) in the uncertainty of that parameter. Corre- 
lations among systematic uncertainties between channels 
are taken into account by using a single parameter for 
the same source of uncertainty. 

We distinguish between systematic uncertainties that 
only affect the yield of signal or background, and those 
that change the shape of the R distribution. We consider 
the jet energy scale, jet energy resolution, jet identifica- 
tion, PDFs, background modeling, and the choice of mt 
in the calculation of P sgn as uncertainties affecting the 
shape of R. Systematic uncertainties on normalizations 
include lepton identification, trigger requirements, uncer- 
tainties on the normalization of background, the uncer- 
tainty on the luminosity, MC modeling, and the determi- 
nation of instrumental background. We also include an 
uncertainty on the templates because of limited statistics 
in the MC samples. 

The statistical and systematic uncertainties on / mea s 
are given in Table [5] We evaluate the size of the indi- 
vidual sources of systematic uncertainty by calculating 
/mcas and a t i using the parameters vy. shifted by ±1SD 
from their fitted mean. 

To estimate the expected uncertainty on the result, 
ensembles of MC experiments are generated for differ- 
ent values of /, and the maximum likelihood fit is re- 
peated, yielding a distribution of / m eas for each generated 
/. Systematic uncertainties are included in this proce- 
dure, taking correlations between channels into account. 
We then apply the "ordering principle" for ratios of like- 
lihoods [28| to the distributions of / m cas and generated 
/, without constraining / mC as to physical values. The re- 
sulting allowed regions for different confidence levels as 



TABLE 2: Summary of uncertainties on /, 



Source +1SD -1SD 



Muon identification 





01 


-0 


01 


Electron identification and smearing 





02 


-0 


02 


PDF 





06 


-0 


05 


mt 





04 


-0 


06 


Triggers 





02 


-0 


02 


Opposite charge selection 


n 
u 


n i 


-u 


fi 1 

Ul 


Jet energy scale 





01 


-0 


04 


Jet reconstruction and identification 





02 


-0 


06 


Background normalization 





07 


-0 


08 


MC statistics 





03 


-0 


03 


Instrumental background 





01 


-0 


01 


Integrated luminosity 





04 


-0 


04 


Other 





02 


-0 


02 


MC statistics for template fits 





10 


-0 


10 


Total systematic uncertainty 





15 


-0 


18 


Statistical uncertainty 





33 


-0 


35 




FIG. 3: (Color online) For all channels the 68.0% (inner), 
95.0% (central), and 99.7% (outer) C.L. bands of / as a func- 
tion of /meas from likelihood fits to MC events. The thin yel- 
low line indicates the most probable value of / as a function of 
/moas, and therefore represents the calibration of the method. 
The vertical dashed black line indicates the measured value 
/moas = 0.74. 

a function of / mC as and / are shown in Fig. [3] From the 
maximum likelihood fit to data, we obtain 

/ moas = 0.74t°;« (stat+syst) . (5) 

The simultaneously extracted ti cross section is found to 
be 

o t - t = 8.3±o'.g (stat+syst) pb (6) 
for mt = 172.5 GeV and in good agreement with the SM 
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prediction of a t i = 1 A%_^^ pb [29(- The comparison of 
/ for prediction and data with the fitted result is shown in 
Fig. [2] The measured fraction is consistent with the SM 
expectation (/ = 1) and we exclude the no-correlation 
hypothesis (/ = 0) at the 97.7% C.L. For the SM value 
of / = 1 we expect to exclude the hypothesis / = with 
99.6% C.L. 

Assuming / moas and using the full matrix elements 
for tt production with SM spin correlation or without 
spin correlation, other observables can be extracted to 
study the impact of this measurement. For illustra- 
tion, we derive C from the measured value of / and 
the NLO prediction of C in the SM, yielding C mcas = 
0.57 ±0.31 (stat+syst) 

In summary, we have presented the first measurement 
of the fraction of tt events with correlated spins us- 
ing a matrix element technique. This fraction can be 



translated into the most precise value of the correlation 
strength C meas to date. 
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